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Abstract. The density function for the joint distribution of the first and second 
eigenvalues at the soft edge of unitary ensembles is found in terms of a Painleve II 
transcendent and its associated isomonodromic system. As a corollary, the density 
function for the spacing between these two eigenvalues is similarly character- 
ized.The particular solution of Painleve II that arises is a double shifted Backlimd 
transformation of the Hasting-McLeod solution, which applies in the case of the 
distribution of the largest eigenvalue at the soft edge. Our deductions are made 
by employing the hard-to-soft edge transitions to existing results for the joint 
distribution of the first and second eigenvalue at the hard edge [13]. In addition 
recursions under a a + 1 of quantities specifying the latter are obtained. A 
Fredholm determinant type characterisation is used to provide accurate numerics 
for the distribution of the spacing between the two largest eigenvalues. 



1. Introduction 



Fundamental to random matrix theory and its applications is the soft edge 
scaling limit of unitary invariant ensembles. As a concrete example, consider 
the Gaussian unitary ensemble, specified by the measure on complex Hermitian 
matrices H proportional to exp(— TrH^)(dH). This measure is unchanged by 
the mapping H 1— > UHU^, for U unitary, and is thus a unitary invariant. To 
leading order the support of the spectrum is {^2N,\/2N), although there is a 
nonzero probability of eigenvalues in (—00, -V2N) U (/2N, 00), and for this reason 
the neighbourhood of\/2N (or -\/2N) is referred to as the soft edge. Moreover, 
upon the scaling of the eigenvalues 1— 7>v^2N + X(A/2N^^^, the mean spacing 
between eigenvalues in the neighbourhood of the largest eigenvalue is of order 
unity. Taking the N — )• 00 limit with this scaling gives a well-defined statistical 
mechanical state, which is an example of a determinantal point process, and 
defined in terms of its A:-point correlation functions by 



det 



K 



soft / 



j,e=i,...,k 



(1.1) 
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where - referred to as the correlation kernel - is given in terms of Airy 
fiinctions by 

r°f'(x,y) := Ai(x)Ai^(y)-Ai(y)Ai^(x) ^ ^^ ^^ 

The determinantal form (i.i) implies that in the soft edge scaled state, the 
probability of there being no eigenvalues in the interval (s, oo), is given by [9] 

fo / ■\\k /.oo /-oo 

E|°ft(0; (s,cx5)) = 1 + E ^-jr- I ■ ■ ■ i dx/c X,), 

= det(l-K-^)), (1.3) 

where is the integral operator on (s,oo) with kernel K^°^^{x,y) (as given 

in (1.2)). The first equality in (1.3) is generally true for a one-dimensional point 
process, while the second equality follows from the Fredholm theory [31]. The 
structure of the kernel (1.2) makes it of a class referred to as integrable [16], and 
generally this class of integrable kernels have intimate connections to integrable 
systems. Indeed one has that [29] 

det (1 - K-f^)) = exp (^- ^"(f - s)ciHt)dt^ , (1.4) 

where q{t) satisfies the particular Painleve 11 ordinary differential equation ( ' = 
d/dt) 

q = 2q^ + tq, (1.5) 

subject to the boundary condition 

q{i)r ^HO- (1-6) 

Our interest in this paper is in the joint distribution of the largest and second 
largest eigenvalue at the soft edge, and the corresponding distribution of the 
spacing between them. Let p^2f (^i'^2)/ ^1 > ^2, denote the density function of 
the joint distribution. Then analogous to the first equality in (1.3) we have 



soft 

\2) 



(xi,X2) = det 



X det 



'K^^\xi,Xi) K^°^\xi,X2) 



E 

k=l 



M 



K^°^\x2,Xi) 



i=i,-,k 



dyi 

X^°f*(xi,X2) 

K'°^\X2,X2) 

K-ft(y^,X2) 



12 



dyic 



j=i,-,k 



K»ft(xi,y,) 
K«°ft(x2,y,) 



f=l,...,k 

e=i,...,k 

j,l=l,...,k. 



(1-7) 
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With A^°^'(s) denoting the density function for the spacing between the two largest 
eigenvalues we have 



isoft/gx 



/oo 
dxp^i°f{x + s,x). (1.8) 



We seek to characterize (1.7) and (1.8) in a form analogous to (1.4). This involves 
functions which are components of a solution of a particular isomonodromic 
problem relating to the PIT equation. Such characterizations have appeared in 
other problems in random matrix theory and related growth processes [1], [13], 
[5], [27]- 

The starting point for us is our earlier study [13] specifying the joint distribu- 
tion of the first and second smallest eigenvalues, and the corresponding spacing 
distribution between these eigenvalues, at the hard edge of unitary ensembles. In 
random matrix theory the latter applies when the eigenvalue density is strictly 
zero on one side of its support, and is specified by the determinantal point process 
with correlation kernel 

[x,y)~ 2{x-y) ' ^^^^ 

where x, y > 0. Note the dependence on the parameter a (a > — 1) which physically 
represents a repulsion from the origin. The relevance to the study of the soft edge 
is that upon the scaling 

X _ 22/3^-2/3^]^ (1.10) 

(and similarly y), as fl ^> 00 the hard edge kernel (1.9) limits to the soft edge kernel, 
and consequently the hard edge state as defined by its correlation functions limits 
to the soft edge state [4]. Thus our task is to compute this limit in the expressions 
from [13]. Moreover, recurrences under the mapping of the latter will be specified. 

In Section 2 the evaluation of the joint distribution of the first and second 
eigenvalue at the hard edge from [13] is revised. This involves quantities relating 
to the Hamiltonian formulation of the Painleve III' equation, and to an isomon- 
odromic problem for the generic Painleve III' equation. Details of these aspects 
are discussed in separate subsections, with special emphasis placed on the trans- 
formation of the relevant quantities under the mapping a + 1. Second order 
recurrences are obtained. In Subsection 2.4 initial conditions for these recurrences 
are specified. Section 3 is devoted to the computation of the hard-to-soft edge 
scaling of the quantities occurring in the evaluation of the joint distribution of 
the first and second eigenvalue at the hard edge. This allows us to evaluate the 
joint distribution of the first and second eigenvalues at the soft edge in terms of a 
Painleve II transcendent and its associated isomonodromic system. 
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2. Hard edge a > joint distribution of the first and second eigenvalues 

2.1. The result from [13]. Let p^2)'^'"{x-[,X2), X2 > x\ denote the joint distribution 
of the smallest and second smallest eigenvalues at the hard edge with imitary 
symmetry. It was derived in [13] that 

„hard, a i „ „ „\ 
P(2) (S-Z,S) 

zV(s -z)''e-'^/4 ( Ar 



42«+3r(fl + i)r(fl + 2)r2(fl + 3) 



exp y [v[r) + 2C(r)]^ {udzV - vdzu) . 



(2.1) 



Here v(s) is the solution of the second-order, second-degree ODE (' = d/ ds) - a 
variant of the tr-form of the third Painleve equation, [13, Eq. (5.25)] 

s^{v")^ -{a + 2)^{v'f + v'{Av' - l){sv' -v) + \a{a + 2)v' - ^a^ = 0. (2.2) 

Important to our subsequent workings is the fact that p^^^'^'''{s) - the probability 
density function for the smallest eigenvalue at the hard edge of an ensemble with 
imitary symmetry - can be expressed in terms of v{s) by [12], [10, Eq. (8.93)] 

= 2^.«r(„-:'i)r(» + 2) (I (•'(" - 1) t) ■ 

To define C(s), introduce the auxiliary quantity = fi(s) according to [13, 
Eq- (5-25)], 

fi + s = 4sv'. (2.4) 

Then, according to [13, Eq. (5.20)], C is specified by 

u'-2 

2C + fl + 3 = si- . (2.5) 

F 

These quantities are closely related to the Hamiltonian variables of Okamoto's 
theory for Pill', as will be seen subsequently. 

The variables m(z;s) and c(z;s) are the components of a solution to the as- 
sociated isomonodromic problem for the generic third Painleve equation or the 
degenerate fifth Painleve equation. They satisfy the Lax pair [13, Eqs (5-34-7)]/ on 
that domain s > z, s, z e R, with real a > — 1, a e M, 

z{s — z)dzU = —Czu — {}i + z)v, (2.6) 

z(s - z)dzv = -z + I (z - s)] M + [-2s + {C + a + 2)z]v, (2.7) 



and 



{s — z)sdsU = zCu + {}i + s)v, (2.8) 
(s — z)sdsV = z£,u — [s(2C + fl) — zC]v, (2.9) 
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where ^ is a further auxiliary quantity specified by ([13, Eq. (5.19)]) 

sC(C + a) 



(2.10) 



}i + s 

For (2.6)-(2.9) to specify a unique solution appropriate boundary conditions 
must be specified. Their explicit form can be found in [13]. 

2.2. Okamoto PHI' theory. We seek to make the links to the Hamiltonian theory 
of the third Painleve equation in order to draw upon the results of Okamoto 
[24], [23] and the work by Forrester and Witte [12]. As given in these works the 
Hamiltonian theory of Painleve 111' can formulated in the variables {q,p;s,H} 
where the Hamiltonian itself is given by ( ' = d/ds) 

sH = q^p^ - {q'^ + viq - s)p + l{vi +V2)q. (2.11) 

With H so specified the corresponding Hamilton equations of motion are 

sq' = 2q^p - {q^ + viq - s), (2.12) 

sp' = -Iqp^ + {2q + vi)p - \ {vi+V2). (2.13) 

From these works its known that the canonical variables can be foimd from 
the time evolution of the Hamiltonian itself by 

p = h' + \, (2.14) 

_ sh" - vih' + \v2 



where 



(2.15) 



h = sH+\v\-\s. (2.16) 
In turn the Painleve 111' u-function is related to the Hamiltonian by 

c^iii'(s) := -{sli)\^^,,^-\vi{vi-V2) + \s. (2.17) 

In the work [13] (see Prop. 5.21) the identification made with the Painleve 111' 
system gave the parameter correspondence Ci — a + 2,V2 = ii — 2 and 

v{s) = -aiii,{s) + \s - a -2. (2.18) 

The quantity C appearing in (2.1) and the auxiliary quantities ji and ^ can be 
related to p and q in the corresponding Hamiltonian system. 

Proposition 1. The variables }i, C, ^ are related to the canonical Painleve HI' co-ordinates 
by 

f = (p-l)L^,/4' (2-19) 
C = -qp\,^,/i, (2.20) 
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Proof. From Eq. (5.21) of [13] and (2.18) we compute that 

i/(s)=/!(s/4)-|(fl + 2)2+is. (2.22) 

Differentiating this and employing the relations (2.4) and (2.14) we find (2.19). 
Using (2.5) we note that As^v" = 2s + (2C + fl + 2)j,i and with the above equation 
and (2.15) we deduce (2.20). Equation (2.21) then follows from (2.10). □ 

For the Hamiltonian (2.11), Okamoto [24] has identified two Schlesinger trans- 
formations with the property 

Ti{vi,V2) = ivi + l,V2 + l), T2ivi,V2) = {vi+l,V2-l), (2.23) 

and has furthermore specified the corresponding mapping of p and q. Recalling 
{vi,V2) in terms of a above (2.18), we see that in the present case Ti corresponds 
to fl fl + 1. Reading from [12] Eqs. (4.40-3) gives the following result. 

Proposition 2 ([12, Eqs. (4.40-3)]). The Painleve IIT canonical variables q[a\ {s),p[a\ (s) 
satisfy coupled recurrence relations in a 

+ ^] = + q[a] {q[a] l^]-\)-2)+s' ^"■"4) 

p[a + 1] = -^q[a] {q[a] {p[a] - 1) - 2) + 1. (2.25) 

The reader should note that we haven't made the scale change s 1— )• s/ 4 here. The initial 
conditions are given by (2.38) below for the sequence a e Z^o- 

2.3. Isomonodromic system. We now turn our attention to the isomonodromic 
system (2.6) -(2.9) for u,v associated with the Painleve system. Following the 
development of [13] we define the matrix variable 



Y(z;s) = rr';M. (2.26) 




To begin with our interest is in the recurrence relations that are satisfied by u 
and V upon the mapping a 1— > a + 1. 

Proposition 3. The isomonodromic components u, v satisfy linear coupled recurrence 
relations in a 

+ 1] = ^ ("W - ^jflP^'I^] j ' (^-^7) 

r 11 s lC[a\+af C[a\+a .A , „s 

v[a + l] = - — 4^ (zu[a]-s^v[a]y (2.28) 

The initial conditions are given by {2.^g)for the sequence a e Z^q- 

Proof. The result (2.1) from [13] was derived as the hard edge scaling limit of 
the joint distribution of the first and second eigenvalues in the finite N Laguerre 
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unitary ensemble. In the latter, the transformation a i— > a + 1 implies a Christoffel- 
Uvarov transformation of the weight zv{x) i— > (x + t)zv{x). From the work of 
Uvarov [30] we deduce that the orthogonal polynomials piq{x; t; a) (we adopt the 
conventions and notations of Section 2 in [13], which should not be confused with 
their subsequent use in Section 3) transform 



pN ■= PN{x;t;a + l) 



,{i,o) 



[pN+l{x;t;a)pN{-t;t;a) - pj^{x;t;a)p]^^i{-t;t;a)], 

(2.29) 

where A^'*^' is a normalisation. In the notations of [13] the three term recurrence 
coefficients transform as 



PN+i{-i; f;fl)pN-i(-f; t;a) 



7n+i7n-i 



PN{-t; t;aY 



bN = bN + ciN 



PN-i{-t; t;a) 
PNi-t; t;a) 



■ «N+1 



PN+i{-i; t;a)' 



(2.30) 
(2.31) 



Employing the variables Qn/^N (see the definitions Eqs. (3.41) and (3.52) of [13]) 
instead of pjv, Pn-i we find that the transformation gives 



Qn{x) 
Rn{x) 



t 

X + t 

t 

X + t 



'<N+l[-t) 



QN{-t) 



Rn{x) Rn+i{x) 



RN{-t) RN+ii-t) 



(2.32) 

(2-33) 



However the second of these equations will suffer a severe cancellation imder the 
hard edge scaling limit t ^ s/ AN, x — z/ 4N as N 00 so we need to be able 
to handle the subtle cancellations occurring. For this we employ a restatement of 
the identity Eq. (3.42) of [13] 

xOnQn + {kn - t)RN - {kn+1 + i)RN+i = 0, (2.34) 
which gives us an exact relation between and Kjv+i- now compute 
t 0nQn(-O tQM{-t)RNix)+xRN{-t)QN{x) 



R^ix) 



x + t RN{-t) - t)RN{-t) - WNQN{-t) 



(2-35) 



We are now in a position to take the hard edge scaling limits, as given in [13] by 
Eq. (5.12) for kn, Eq. (5.10) for 9}^, Eq. (5.28) for Qjv and Eq. (5.29) for _Rjv. In 
addition we employ the identity, [13, Eq. (5.45)] 



v{s;s) _ g _ sC 
u{s;s) C + a ji + t 



(2.36) 



The final result is (2.27,2.28) where all dependencies other those other than a are 
suppressed. □ 



8 



N. S. WITTE, F. BORNEMANN, AND P. J. FORRESTER 



2.4. Special Case a E Z. Jn Subsection 5.2 of [13] determinantal evaluations of 
the Painleve variables v, }i, C and ^, of the isomonodromic components u, v and 
Aa for a E Z^o were given. These were of Toeplitz or bordered Toeplitz form 
and of sizes a x a, {a + 1) x {a + 1) and (a + 2) x (a + 2) respectively. Here we 
content ourselves with displaying the first two cases only, which can serve as 
initial conditions for the recurrences in Propositions 2 and 3. In order to signify 
the fl-value we append a subscript to the variables. In all that follows Ia-{z) refers 
to the standard modified Bessel fimction with index a and argument z, see §10.25 
of [25]. 



2.4.1. = 0. Some details of the first case a = were given in Propositions 5.9, 
5.10 and 5.11 of [13] and we augment that collection by computing the remaining 
variables. Thus we find for the primary variables 

vo{s) = 0, /io(s) = -s, Co(s) = 0, ^o(s) = 0, (2.37) 

for the canonical Hamiltonian variables 

po(s/4)=0, ^o(s/4)=^|^, (2.38) 

the isomonodromic components 

8 4 

uo{z;s) = -hi^/z), vo{z;s) = (2.39) 
z ^yz 



and the distribution of the spacing 



1 



-z/4 



(2.40) 



This formula is essentially the same as the gap probability at the hard edge 
for fl = 2, as one can see from the ^ = specialisation of Eq. (8.97) in [10]. 
Interestingly we should point out that the moments of the above distribution can 
be exactly evaluated and we illustrate this observation by giving the first few 
examples {mo = 1) 

mi =4e2 [Io(2)-Ii(2)], 

m2 = 32e2jo(2), 

m3 = 384e2 [2Io(2) + Ii(2)] , 

m4 = 2048e2[l3Jo(2)+9Ji(2)], 

OT5 = 20480e2 [55Io(2) + 42Ii(2)] , 

me = 98304e2 [557Jo(2) +441Ii(2)] . 
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2.4.2. a = 1. This case was not considered in [13]. We have computed these 
from the results for the finite rank deformed Laguerre ensemble, as given in 
Section 4 of [13], and then applied the hard edge scaling limits given by the Hilb 
iype asymptotic formula Eq. (5.2) therein and the limits of Proposition 5.1 and 
Corollary 5.2 of [13]. For the primary variables we find 

4 4/3 (^/s)2 2/3 (/s) 



the Pill' canonical Hamiltonian variables 



p,(s/4) = l-W^, (..43) 



''^'^'^ = 2WT) — wmM^hm — ' ^"''^ 

the isomonodromic components for generic argument s > z > 

ui[z,s) = — , (2.47) 

Z-fsVS) S-2 

Vsl2(/S)7sl2^)i3^) -Vzh{/Z)hi^s) 
V\{Z,S} — , I2.4»j 

^Jzlo,ysY S-Z 
and the isomonodromic components on s = z 

Ml (s; s) = + r-.\ (2.49) 



and the distribution of the eigenvalue gap is 



Ai(z) = 2-4 [h(^z)h^z) - h{^zf\ ds e-''%irs) 



f<Xi 



s — z 

(2.51) 



From the point of view of checking one can verify that the above solutions satisfy 
their respective characterising equations. 
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2.5. Lax Pairs. We now examine the isomonodromic system from the viewpoint 
of its characterisation as the solution to the partial differential systems with respect 
to z and s. 

Proposition 4 ([13, Eqs. (5.51,5.52,5.54-7)]). The matrix form of the spectral derivatives 
(2.6) and (2.7) and deformation derivatives (2.8) and (2.9) yield the Lax pair 



0^ 

1 0, 



and 



-C 



}i + s \ 
-C — a I 



-C, 



}i + s \ 
-C — a , 



(2.52) 



(2-53) 



This system is essentially equivalent to the isomonodromic system of the fifth Painleve 
equation but is the degenerate case. The system has two regular singularities at z = 0,s 
and an irregular one at z = 00 with a Poincare index of j. 

The form of the isomonodromic system (2.7) -(2.10) is not suitable for comput- 
ing the hard-to-soft edge scaling limit, so we need to perform some preliminary 
transformations on it. 



Proposition 5. Under the gauge transformation 

u,v^ z-^s'''^{s-zy'^u,v 
the spectral derivatives (2.6) and{2.'y) become 



z(s — z)dzU 
z(s — z)3zC 



-z{C+\a) u—{pi + z)v, 



(C- 



\a)z 



(2.54) 

(2.55) 
(2.56) 



whilst the deformation derivatives (2.8) and (2.9) become 

(s — z)sdsU = (C + 2^)zu + {pi + s)v, (2.57) 

(s - z)sdsV = z^u + (C + 2fl)(z - 2s)v. (2.58) 

Furthermore let us scale the spectral variable z sr. Consequently Equations (2.55) and 
(2.56) become 







C + a/2 
[ ? -C-a/2^ 

and Equations (2.57) and (2.58) become 



-C-fl/2 



^ 
-C-a/2, 



C + a/2 



Y, (2.59) 



T. (2.60) 



This system has two regular singularities r = 0,1 and an irregular one at r = 00 with 
Poincare rank of j (due to the nilpotent character of the leading matrix in 2.59), and is 
denoted by the symbol (l)^^)- 
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The precise solutions we seek are defined by their local expansions about 
r = 0,1, and in particular the former case. From the general theory of linear 
ordinary differential equations [28], [6] we can deduce the existence of convergent 
expansions about r = (s 7^ 0) 

00 

u{r;s) = J2 u„,(s)r^o+"', (2.61) 

m— 
00 

i;(r;s)= J]v,„(s)r^o+'«, (2.62) 

m=0 

with a radius of convergence of at most unity. The indicial values Xo are fixed by 

f/vo + suo (a;o - 1 ) = 0, vo (;3Co + 1 ) = 0, (2.63) 

and the appropriate solution has vq = 0, ug 7^ and Xo = 1 (actually from [13, 
Eqs. (5.38, 5.39)] we know uq = s, vq = 0). The general coefficients are given by 
the recurrence relations 

4?n(7fi + 2)sum = - [-4:{m^ + m- 2)s + 2(ot + 2)s(2C + a) + }i{s - 4^)]u,„_i 

+ s/iu„_2 - 2[2(m + 2)s + (2C + a)ii + 2{m + l)/i]v,„_i, (2.64) 



4(ot + 2)v„, = (s - 4^)u„,_i - su„_2 + 2[2C + a + 2{m + l)]v„_i. (2.65) 
The first few terms are given by 

Uo = s, vo = 0, (2.66) 

ui = -i(2C + fl)s + i^((^-|s), VI = -is(?- is). (2.67) 



Similar considerations apply to the local expansions about r = 1 however in the 
hard-to-soft edge limit this singularity will diverge to 00 and we will not be able 
to draw any simple conclusions in this case. 

3. Hard to Soft Edge Scaling 

The hard edge to soft edge scaling limit [4] will be interpreted as the degenera- 
tion of PUT to Pll. Therefore we begin with a svimmary of the relevant Okamoto 
theory for Pll. 

3.1. Okamoto PII theory. Henceforth the canonical variables of the Hamiltonian 
system for Pll will be denoted by {q,p;t,H} and should not be confused with 
the use of the same symbols for PUT. Conforming to common usage we have the 
parameter relations a = a-[ — ^ = 2 — ^o- The Pll Hamiltonian is 

H = -lilq^ - p + t)p - a^q , (3.1) 

and therefore the Pll Hamilton equations of motion (' = d/dt) are 

q = p-q^-\t, p = 2qp + Cii. (3.2) 
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The transcendent Cj{t;(x.) then satisfies the standard form of the second Painleve 
equation 

cj = 2q^ + tq + ci. (3.3) 

The PIT Hamiltonian H(f) satisfies the second-order second-degree differential 
equation of Jimbo-Miwa-Okamoto a form for Pll, 

{Hf+A{Hf+2H[tH-H]-lal = 0. (3.4) 

Using the first two derivatives of the non-autonomous Hamiltonian H 

n = -\p, 

H = -qp - jtti , 
we can recover the canonical variables of (3.2). 

3.2. Degeneration from PHI' to PII. We know from [4] that upon the scaling 
(1.10) of the variables and taking a ^ co the hard edge kernel (1.9) limits to 
the soft edge kernel (1.2) and furthermore the joint distribution p^^'^''^{xi,X2) 

limits to ppj*(xi,X2)- The same holds true for the relationship between p^^^'^''^[s) 
and p^°|*(s). This latter fact helps in our computation of p^°|'(xi,X2)/ since the 
evaluation of p|°j'(s) in terms of Pll is known from previous work [10], [11] 
allowing the limiting form of v[t) in (2.3) to be deduced. But v{t) is the very same 
PUT quantity appearing in the evaluation (2.1) of p^'^'''{s — z,s). 

Proposition 6. Let 

S=fl2[l_22/V2/3T]. (3.6) 

We have that for a — )• 00 

" I + ^ ^ ~ (I) (3-7) 



where t = —2^^^t, 
and furthermore 



C7„(T) = -2l/3H(f)L^=2, (3.8) 



'^n(T)^:^^^logX-^T,T). (3.9) 



Proof. We know from [10, Eq. (8.84)] that 

P(?f (s) = Paf(s) exp (^- (^auit) - ^ logp-f (O) dt^ , (3.10) 

where |0^°|*(s) = _K^°^(s,s). On the other hand 

(s) = lim 22/3fl4/3 hard,«(^2H _ 22/3^-2/3gn^ ^^^^^ 
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Substituting (2.3) in the RHS and (3.10) in the LHS of (3.11), and comparing the 
respective large s forms implies (3.7). The boundary condition (3.9) is immediate 
from (3.10). □ 

It will be shown in the Appendix that the solution of the cr form of PII (3.4) 
with ai = 2 as required by (3.8), and subject to the boundary condition (3.9), can 
be generated from the well known Hasting-McLeod solution of PII. 

The scaled form of C in (2.1), as well as the auxiliary quantities pi (2.4) and ^ 
(2.10) can now be found as a consequence of (3.7). 

Proposition 7. Let s be related to t by (3.6), and define t = —l^^^r as before. As a ^ 00 
F(s) ^ -2i/V/3p(i), (3.12) 

2C(s)+fl^22/V/3 



f(s)^l«2_ 2-2/3^4/3 



^(0 + -77Y -2PW 



(3-13) 
(3-14) 



Proof. Simple calculations using (3.7) and (2.4), (2.5) and (2.10) give (3.12), (3.13) 
and (3.14) respectively. □ 

Now we turn to task of deducing the appropriate scaling of the associated 
linear systems and their limits in the hard edge to soft edge transition. There are a 
handful of references treating the problem of how the degeneration scheme of the 
Painleve equations is manifested from the viewpoint of isomonodromic deforma- 
tions. In comparison to the work [18] our situation is that of the degenerate PV 
case with nilpotent matrix Aco as given by Eq. (11) in that work and its reduction 
to the case of Eq. (13), again with nilpotent matrix Aoo, which corresponds to PII. 
In the more complete examination of the coalescence scheme, as given in [21], our 
reduction is the limit of the degenerate PV (P5-B case) to that P34, and therefore 
equivalent to PII. However many details we require are missing or incomplete in 
[18] and [21], so we give a fuller accoimt of this scaling and limit for our example. 

Lemma 1. Let e = 2^^^a~^^^. The independent spectral variable scales as r = e^x. 
Under the hard-to-soft edge scaling limit e the isomonodromic components scale as 
u = 0{e-^) and v = 0{e~^). 

Proof. Let us denote the leading order scaling of the expansion coefficients given 
in (2.61), (2.62) by u,n = Ola^"') and v,„ = 0{a'^"'). Employing the leading order 
terms of the auxiliary variables (3.12), (3.13) and (3.14) in the recurrence relations 
for the coefficients (2.64) and (2.65) we deduce that 

COm = max{a;„,_i + ^,CVm-2 + J/^m-l}/ 

Am = max{a;m_i + i,cOm-2 + 2, A^-i + i}- 
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In fact all terms on the right-hand side balance each other and are satisfied by the 



single relation cOm = co 



m—l 



Am-1- The solution to these is co„ 



■2, 



A 



jOT + J, given the initial condition oJi = |. We then deduce that each 
term in the expansions has leading order u„,r"'+^ = 0{e^^), Vm''"'^^ = 0(e^^), 
independent of m. Given that the expansions converge uniformly then the whole 
sums have the stated leading order expansions. □ 

Proposition 8. Let the isomonodromic components scale as u{r;s) = U{x;t), v{r;s) = 
e~-^y (x; f), as only the relative leading orders matter. As e — >■ the spectral derivative 
scales to one of the Lax pair for the second Painleve equation t,x ^M. 





-I 0, 



X + 



l{t-p) + 




and the deformation derivative scales to 





(3-15) 



(3.16) 



Proof. Our starting point is the Lax pair for the degenerate fifth Painleve system 
given in Equations (2.59) and (2.60). Using the expansions (3.13,3.12,3.14) and (3.6) 
we deduce that the matrix elements appearing in this pair scale as 

C + fl/2 



-C-a/2 + 



1 



-e ^ho, 



-C - a/2 — ~ ihox, 

v I 



1 C + fl/2 _2 

- H ^ ~ e 

r r — 1 



jl + S 



(C + fl/2) 



2-ifl2 



1 



r-1 
}i + s 1 



\{t-x) + \hl-\^ 



1, 



s r — 1 

sr s 
where the abbreviation is 



17) 
18) 
19) 
20) 
21) 
22) 



hQ=2\q + 



(3-23) 



Using the scaling for u, v we deduce a meaningful limit as e ^> given in Equa- 
tions (3.15) and (3.16). One can check that the compatibility of these two equations 
is ensured by the requirement that q, p satisfy the Hamiltonian equations of motion 
(3.2). □ 
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Remark i. For general ai the Lax pair of the PIT system is 



Iit-P) + 





and 



dtY 





(3-25) 



Equation (3.24) has the same form as the nilpotent case of Eq. (13) of [18], and in 
addition both members of the Lax pair (3.24,3.25) are a variant of the system Eq. 
(31), given subsequently in [18]. It has also been shown in [19] that the Lax pair of 
this latter system is related to that of Flashka and Newell [7] via an "unfolding" of 
the spectral variable supplemented by a gauge transformation (see also 5.0.54,5 
on pg. 175 of [8]). In contrast the Lax pair of Jimbo, Miwa and Ueno [17] is 
not equivalent to any of those mentioned above. In the hard-to-soft edge scaling 
the regular singularity at r = has transformed into the regular singularity at 
X = 0; the regular singularity at r = 1 has merged with the irregular one at r = 00 
yielding an irregular singularity at x = 00 with its Porncare rank increased by 
unity, now being |. The symbol of the new system is (1)(2). 



The solution we seek can be characterised in a precise way though its expansion 
about the regular singularity x = 0. 

Lemma 2. Let us assume \ p{t)\ > S > and t,q{t),p{t) lie in compact subsets ofC. 
The isomonodromic components U, V have a convergent expansion about x = 0, with 
indicial exponent xo = 1/ whose leading terms are 

p^{2q^ - p + t) +2qp - A 2 



U{x;t) = 2x ■ 



3p 

p {p^{2q^ -p + t)+Apq-8) {2q^ - p + t) + 2{p^ - 16q - Apt) 



A8p 



V{x;t) 



+ 0{x^), 
p^{2q^-p + t)+ 8pq + 8 ^2 



(3.26) 



+ 



3p2 

p {p^{2q^ -p + t)+ \2pq + 24) {2q'^ -p + t)+ 2{5p^ + 16q 



+ 0(x4). 



(3-27) 



Proof. The scaling relations (3.6, 3.7, 3.12, 3.13, 3.14) can be applied term-wise to 
the expansions (2.61, 2.62) along with the explicit results for the leading coefficients 
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(2.66, 2.67). Alternatively one can compute the recurrence relations for the local 
expansion of the system (3.15) 

00 

U,V= Y^\Jm,Y„xX0+'". (3.28) 

m=0 

In such an analysis one finds for the leading relation p [{xo ~ l)Uo — pVp) = 
and 2p^{xo + l)Vo = 0. Clearly for a well-defined solution at x = we must have 
Vq = and so xo = ^ with Uq 7^ 0. The recurrence relations for the coefficients 
are given by 

2m{m + 2)pUm = [p^{2q^ -p + t) + {A- 2m)qp - 4ot]U„_i 

+ 2p{pq - m)V„_i - p^U,„_2, (3-29) 



2{m + 2)p2v,„ = [p2(2^2 _ p ^ + + S\\J„,-i 

+ 2p{pq + 2)V„,_i - p2u„,_2. (3-30) 

Given Uq = 2 these recurrences generate the unique solution stated above. From 
these recurrence relations it is easy to establish that the local expansions (3.28) 
define an entire function of x. □ 



We now have all the preliminary results to obtain the sought Painleve II 
evaluation of p^j'/ specified originally as the Fredholm minor (1.7). 

Proposition 9. Let p^°j'(0 be given by (3.10). For some constant Cq still to be determined, 
and boundary conditions on U and V still to be determined 

pl-f{t,t-x) = 

Cop-(0 .-/^exp exp (/;^^dy { (2,+ ^) (-y) "^2^ - A}) 

X (!J3.,.y - Vd:cU) {-2}l^x; -2}'^t). (3.31) 

Frooj. Applying the gauge transformation (2.54) to (2.1) and absorbing the pre- 
f actors exp(— s/4) and s" into the integral we have 

p^^^^-\s - z,s) = a (so) exp ( ^ \v{w) -'^+a + 2C{w)+a 



So 



w 



4 

X (udzV — vdzu) (z;s). 



where Cfl(so) is a normalisation independent of s,z but dependent on a and the 
reference point sq. We are now in a position to apply the limiting forms (3-7,3-i3) 



JOINT DISTRIBUTION OF . . . AT THE SOFT EDGE OF UNITARY ENSEMBLES 17 

to this, thus obtaining 

= Coexp ^|°^^dy (^H + 2q+^^ (-y)^ {Ud,V - VdM) {-2^^^x; -2^/^). 

(3-32) 

To proceed further, we make use of (3.8) and (3.10) to note that for suitable Cq, 
Furthermore, since 



and 



pg'(s) =K-ft(s,s)^^^-i-exp(-4s3/2), (3.35) 



we must have 



y->oo y 

The Hamilton equations (3.5) then imply 



H{-y)^^^J2^+-. (3.36) 



2 2 



'?(-y).,~^-\/?-^' (3-38) 



y^~y y y^ 

'y^oo~\/2 ~ 4y' 
and thus 

5 



2y- 

Consequentially, for suitable Cq, 

Substituting (3.33) and (3.39) in (3.32) gives (3.31). □ 

It remains to specify Cq in (3.31), and furthermore to specify the x,t ^ co 
asymptotic form of U and V. For this we require the fact, which follows from (1.7), 
that for t,x,t — X ^ 00, 

pg'(U-x)~p|?f(Opg'(t-x). (3.40) 
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Proposition lo. In (3.31) 



Co 



Atz' 



and furthermore, for x,t,t — x ^ 00 we have 



with 



U{—x; —t) ^ a{x, t) exp 
V{-x;-t) ^ b{x,t)exp 

f5/4 

a{x, t) 



V2 



i3/2 . 



i3/2 . 



(t-X 



,3/2 



{t-x) 



3/2 



(t-x) 



1/4' 



b{x,t) 



t — X 



a[x, t). 



(341) 

(342) 
(343) 

(344) 
(345) 



Proof. Substituting (3.40) in the LHS of (3.31) and making use of (3.34) and (3-35), 
it follows from (3.31) that in the asymptotic region in question 



Sn{t-x) 



exp 



{-tit--?") 



Cor5/2g^p /_4f3/zj {ud^V -Vd:,U){-2"'x;-2''H). (3.46) 



.4x3/2 



1/3^. _ol/3. 



(347) 



We see immediately from this that (342) and (3.43) are valid, for a(x,t), h{x,t) 
algebraic functions in x and t satisfying 

_L_^ ^ Co \a{2^'\,2^'^t)dM-y,2^'h) 

- b{2^'^x,2^'^t)dya{-y,2^'^t) 
On the other hand, we read off from (3.15) that for x, i ^> 00 

dyU{y;-t)\y_^^ _ U{-x;-t)-V{-x;-t). 

Making use of the leading terms in (3.37) and (3-38) as well as (3.42), (3.43) it 
follows that (3.45) holds true. This latter formula substituted in (3.47) implies, 
upon choosing Co according to (3.41), that 

\Jt — X 

and (3.44) follows upon taking the positive square root. 

We can offer a refinement on the above argument by examining in more detail 
the isomonodromic system in the asymptotic regime t ^ —00. In this regime the 
leading order of the differential equations (3.15) and (3.16) become 



(348) 
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and 



(349) 



In retaining only the leading order terms we have also implicitly assumed that 
X 00 because in the regime as t ^ —00 our solution p is vanishing, which is the 
location of the only non-zero, finite singularity. Clearly 







(3-50) 



so that U{x;t),V{x;t) ~ f{i)gi,2{i — x). This means that dtf = —y—^f with a 
solution proportional to exp ^{—t)^^"^. For the remaining factor we have 



Si 
yS2j 



or, in terms of the components 



gl=lix-t)gi, g2=g'l 



t 



(3-51) 



(3-52) 



Thus, with Ai and Bi the two linearly independent solutions of the Airy equation 
[25, §9.2], we have 

Uix;t) ^ exp ^(-0^/2 {«(0Ai(2-i/3(x - t)) + /5(0Bi(2-i/3(^ _ j)) j , (3.33) 
and 

y(x;0~exp^(-03/' 



-m 



2-1/3^1' (2-i/3(x - f)) + W-iAi(2"i/3^x - 0) 



2-i/3g./^2-i/3(;c-0) + 



Bi{2~^/\x - t)) 



(3-54) 



From these we compute 



UdrV - VdrU ~ 



expl'^i-tf/^ {-lix - t) [aAi + /3Bi]2 +2-^/3 [^Ai' + fiBif} . (3.55) 

Clearly /5 = in order to suppress the dominant terms, and by employing the 
exponential asymptotics of the Airy function through to second order (the leading 
order cancels exactly) we have 



2-4/3^2 

{Ud,V - Vd,U) (x; t) ~ exp 



i-^i-tr^-i-^ix-tf'^ 



■ (3-56) 
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Employing this into the factorisation of p^2f(^' f — ^) as t — x, x — > oo we deduce 

Co«(-2l/3f)2=22/3i5/2, (3.57) 

and with Cq = l/in we infer a{t) = 2^^^^^7t^^^( — f)^'''*. This then precisely 
reproduces the boundary conditions given by (3.42, 3.43) along with (3.44, 3-45)- 
A by-product of this argument is that the isomonodromic components can be 
represented in the forms 

, Ai(2-l/3(;C_f)) 

Aii-2^UH) ' (3-58) 



2-i/3Ai'(2-i/3(x - 0) + A/-^Ai(2-i/3(x - t)) 
V{x; t) ^ t Ai(-2-i/30 ' ^3.59) 

although this is still only valid in the regime t — x, x— )•— 00. □ 



We remark that as written (3.31) is not well defined for f < 0. In this region we 
should use instead (3.32), with a suitable fo to make use of (3-31) for f > so that 
Co can be specified. 

In [13] the analogue of (3.31)/ Eq. (2.1), was used to provide high precision 
numerics for the spacing between the two smallest eigenvalues at the hard edge. 
But to use (3.31) to compute the spacing distribution (1.8) presents additional 
challenges to obtain control on the accuracy. The essential problem faced in 
comparison to [13], and in also in comparison to the study [26] relating to PDEs 
based on the Hasting-McLeod PII transcendent, is that our boundary conditions 
involve algebraic terms, and thus cannot be determined to arbitrary accuracy. In 
relation to p{—y) and ci{—y) this can perhaps be overcome by using the theory of 
the Appendix to map to the Hasting-McLeod solution. But even so, the problem of 
extending the accuracy of the algebraic terms a{x, t) and h{x, t) in (3.42) and (3.43) 
would remain. While these points remain under investigation, the problem of 
determining numerics for (1.8) can be tackled by using a variant of the Fredholm 
type expansion (1.7), as we will now proceed to detail. 



4. Numerical Evaluation of Moments 

We provide some numerical data for the distribution of the spacing of the 
two largest eigenvalues based on the accurate numerical evaluation of operator 
determinants that is surveyed in [2]. The joint probability distribution of the two 
largest eigenvalues is amenable to this method since it is given in terms of a 2 x 2 
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Table i. The first four statistical moments of the density function A'*°''(s). 
mean variance skewness excess kurtosis 



1 .90435 049 0.68325 206 0.56229 2 



0.27009 



operator matrix determinant: 



F{x,y) 



5fft(0;(y,oo)) 



dz 



det 



/ 



V 



2 ^soft j^soft 
2 ^soft j^soft 



(x ^ y), 
{x > y). 



|L2(y,.r)©L2{.T,oo)/ 

Here, the differentiation with respect to z can be accurately computed by the 
Cauchy integral formula in the complex domain [3]. It has been used in [2] to 
calculate the correlation between the two largest eigenvalues to 11 digits accuracy 

p = 0.50564 723159. 

This number can be calculated without any differentiation of the joint distribution 
function F since, according to a lemma of Hoeffding [15], the covariance is given 
by 

/OO rOO 
/ {F{x,y) — F{x,oo)F{oo,y)) dx dy. 
-00 J — CO 

In contrast, the spacing distribution 



G(s 



rS roo J 

(s) = A-f'(c7)d^ = J_^-F{x,y)U=y+sdy 

requires numerical differentiation in the real domain which causes a loss of a 
couple of digits. The differentiation is done by spectral collocation in Chebyshev 
points of the first kind and G is numerically represented by poljmomial inter- 
polation in the same type of points; Table 2 tabulates the values of A^°^'(s) and 
G(s) for s = 0(0.05)8.95 to an absolute accuracy of 8 digits based on a polynomial 
representation of degree 64 that is accurate to about 9 to 10 digits. Fig. 1 plots the 
density function as compared to a histogram obtained from 10 000 draws from a 
1000 X 1000 GUE at the soft edge. 

The moments of the random variable S representing the spacing are obtained 
from the following derivative-free formulae obtained from partial integration: 

/■OO I^CO 

E(S")=/ s"dG(s)=n s"-^(l-G(s))ds (n = l,2,...). 
Jo Jo 

This way we have obtained the first four statistical moments shown in Table 1; 
estimates of the approximation errors by calculations to higher accuracy indicate 
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spacing s 

Figure i. A plot of the density function A^°^*(s) as compared to a 
histogram of 10000 draws from a 1000 x 1000 GUE at the soft edge. 

the given digits to be correctly truncated. The total computing time was 5 hours 
for the solution and 30 hours for the higher accuracy control calculation. 
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Appendix 

As a technical matter we will need to make use of the Gambler or Folding 
transformation for PII. The fundamental domain or Weyl chamber for the PII 
system can be taken as the interval a E ( — 1,0] or a E [0, j), and there exist 
identities relating the transcendents and related quantities at the endpoints of 
these intervals. In particular, denoting the transcendent q{t;a) and with = 1, 
t = —2^^^s we have [14] 



;q{s;0) . 



(A.i) 



q{s;0) ds 

In addition we will employ the Backlimd transformation theory of PII as 
formulated by Noumi and Yamada (see [20]) and put to use in the random matrix 
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Table 2. Values of the probability density and distribution function for 
the spacing s between the two largest eigenvalues; with s — 0(0.05)8.95. 
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context by [ii]. We define a shift operator corresponding to a translation of the 

(1) 

fiindamental weights of the affine Weyl group ', 

T2 : ao 1-^ ao — 1, 1-^ «! + 1 . (A. 2) 

The discrete dynamical system generated by the Backlimd transformations is also 
integrable and can be identified with a discrete Painleve system, discrete dPI. 
The members of the sequence {^[n]}5^^o/ generated by the shift operator T2 with 
the parameters {olq — n,a.i + n), are related by a second-order difference equation 
which is the alternate form of the first discrete Painleve equation, a-dPl, 

a+l+n a— 1+n 

+ 7T— iV^ = -VM-i- (A.3) 



^[n] +(/[« + !] q[n-l]+q[n] 
The full set of forward and backward difference equations are [22] 

ci[n + l] = -q[n]-^^^^ , (A.5) 
p[n-l] = -p[n]+2q[nf + t , (A.6) 

p[n + l]=t-p[n]+2(q[n] + ^ + J + . (A.7) 



In addition one should note that H[n + 1] = H[n] — q[n + 1]. 

Proposition 11. The solution of the second Painleve equation as given by (3.4) with 
parameter =2 and boundary condition (3.36) is generated from the Hastings-McLeod 
solution by application of the T2 Schlesinger transformation applied twice and the Gambler 
transformation (A.i). 

Proof. Firstly we recall that the parameter for the Hastings-McLeod solution is 
a = 0, ai = 1/2 whereas we have the case of a = 3/2, = 2. Let t = —2~^^^t. 
The leading, and defining, asymptotics of the Hastings-McLeod solution at a = 
0, ai = 1/2 as T ^> +00 is (for f = 1, Eq. (9.47) of [10]) 

q{T;ai = 1/2) ~ Ai(T). 

T— 7-00 

Using the inverse Gambler transformation (A.i) with e = — 1 we have the solution 
as 



q{t;oci=0) ~ -2 



_,/, Ai'{-2-y^t) 



Ai(-2-i/3f) ' 
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and therefore p(f;ai = 0) ^ and H{t;a.i = 0) ^ in this regime. Now using the 
Schlesinger transformations (A.5,A.7) we deduce 

1 

H(f;ai = 2) = H{t;ai = l)+q{t;ai = V 
= H{t;oii = 0) 



= 1)' 

1 



2[q{t;cci = 0)]^ - p{t; oci = 0) + t' 



2[<7(t;«i =0)]2 + t' 



Ai(-2-i/3f) 



[Ai'(-2-i/3f)]2^2-i/3f [Ai(-2-i/30]^' 



which is asymptotically equivalent to (3.36). 



□ 
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